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A SHORT PROOF THAT THE NUMBER OF DIVISION STEPS 
IN THE EUCLIDEAN ALGORITHM IS NORMALLY 
DISTRIBUTED 

IAN D. MORRIS 


Abstract. D. Hensley showed in 1994 that the number of steps taken by 
the Euclidean algorithm to find the greatest common divisor of two natural 
numbers less than or equal to n follows a normal distribution in the limit as 
n tends to infinity. V. Baladi and B. Vallee subsequently gave an alternative 
proof for both the classical Euclidean algorithm and several of its close variants, 
based on a detailed investigation of spectral properties of the transfer operator 
associated to the Gauss map, building on deep results of D. Dolgopyat. In this 
article we give a much shorter, albeit less quantitative, proof of this result 
using only basic spectral properties of the transfer operator together with the 
method of moments and a Tauberian theorem due to H. Delange. 

MSC codes: Primary 11A05, 11K50, 68W40; Secondary 37D20, 37D35. 


1. Introduction 

The Euclidean algorithm for determining the greatest common divisor of two 
natural numbers is believed to be the oldest nontrivial algorithm which remains 
in common use to the present day [HI p.335]. By the late 16 th century it had 
been noticed that the slowest case of the algorithm is that in which the input is 
a pair of successive Fibonacci numbers [22]. Rigorous analysis of the number of 
division steps required by the algorithm began in earnest in the 19 th century when 
the first logarithmic upper bound was given by P.-J.-E. Finck (®, for a detailed 
historical overview see [233). Analysis of the number of division steps typically 
required to execute the algorithm did not develop until the second half of the 20 th 
century, when H. Heilbronn m and J. D. Dixon [7] independently showed that the 
average number of divisions required to process a pair of integers both bounded 
by n is asymptotic to ^ log 2 log n in the limit as n —> oo. A further milestone 
was achieved by D. Hensley m who in 1994 showed that the deviation from this 
mean is asymptotically Gaussian for large n, giving also a local limit theorem and 
an upper bound for the rate of convergence to the limit distribution. 

Analyses of variant versions of the Euclidean algorithm were undertaken in the 
late 20 th century using somewhat disparate methods [3HiH6]. B. Vallee [25] sub¬ 
sequently unified the analysis of the average number of division steps of a range 
of Euclidean-type algorithms into a single technique based on the thermodynamic 
formalism of the Gauss-Kuzman-Wirsing operator (also called the density trans¬ 
former or transfer operator) and its close relatives. This approach was extended in 
joint work with V. Baladi (in, see also 01) where Hensley’s result was generalised 
to prove that the number of division steps asymptotically follows a Gaussian distri¬ 
bution throughout a class of variant Euclidean algorithms, with detailed estimates 
on the rate of convergence. Baladi and Vallee’s results also generalised Hensley’s 


l 


2 


IAN D. MORRIS 


result to different cost measures of the algorithm including the binary length of the 
divisor and the number of occurences of particular quotients. 

Both the proof given by Hensley El and that given by Baladi-Vallee [I] are 
long and technically involved. Baladi and Vallee’s treatment in particular rests 
on a delicate investigation of the norm of certain complex perturbations of the 
transfer operator using techniques developed by D. Dolgopyat [3 in the application 
of thermodynamic formalism to geodesic flows. This may lead the reader to wonder 
whether the fact that the number of steps obeys a Gaussian distribution is amenable 
to a shorter proof. In this article we will show that the existence of a Gaussian 
limit distribution for the number of steps (but without information on the rate of 
convergence) can be recovered via a relatively simple strategy using only qualitative 
features of the spectrum of the transfer operator: indeed, the properties which we 
require do not significantly extend beyond those used in Vallee’s earlier article [25] . 
We also hope that the approach presented here may facilitate the generalisation of 
this type of result to related contexts such as umm where the corresponding 
Dolgopyat-type estimate is unavailable or may even be impossible. 


2. Statement of results 

The standard Euclidean algorithm, starting from a pair of integers (u, v) with 
0 < u < v, iterates the following step, which we call the division step: map the 
pair (u,v) to the pair (r,u), where v = qu + r and 0 < r < v. This procedure 
is repeated until a pair of the form (0, k) is obtained, whereupon the integer k is 
returned as the GCD of u and v. The algorithm may be modelled via the following 
parallel procedure: given a rational number x := u/v G (0,1], we map that number 
to the new rational number T(x) := - — |_-J, repeating until zero is obtained. Since 
at each stage the integer r = v — qu equals precisely v — u \y /u \, it is clear that 
both procedures involves the same number of steps when applied to each pair (it, v). 
To measure the number of steps taken to reduce the pair (u, v) is thus equivalent 
to finding the least integer N(u/v) > 1 such that T N ^ u / v \u/v) = 0. As in the 
work of Vallee and Baladi, in addition to studying the integer N(u/v) we are easily 
able to adapt our arguments to the study of more general “cost” measurements, 
which assign to the pair (u, v ) a cost C(u/v ) depending on the particular sequence 
of divisions performed - or, equivalently, on the digits gu > 1 arising in the finite 
continued fraction representation u/v = [a±,... ,a n ). In principle our method is 
capable of studying the case in which the cost c(n ) associated to division by the 
integer n increases logarithmically with n, but in this note we have restricted our 
attention to bounded cost functions c: N —>• [0, +oo) so as to simplify the exposition. 

Let £: [0,1] —> R denote the function £(x) := (i 0<r 2 )(i+:r) w hich is the density 
of the unique absolutely continuous T-invariant probability measure on [0,1]. (In 
subsequent sections we shall also write £ for the extension of that function to a 
certain disc in the complex plane.) For each n > 0 let us define fi n to be the set of 
all pairs of integers (u, v ) G N 2 such that 1 < u < v < n, and let fi n denote the set 
of all (u,v) £ fi n such that gcd(u,u) = 1. Let P n (resp. P n ) denote the uniform 
probability measure on fi n (resp. fi n ). We prove the following result: 

Theorem 1. Let c: N —> [0,+oo) be a bounded function which is not identically 
zero, and for each u/v = [a \,..., a n \ G (0,1) D Q with a n > 2 define C(u/v ) := 
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Sfc=i c ( a k)- Then there exist constants fi, a 1 > 0 such that for all t £ 


lim P„ 

n—>oo 


(u,v) £ 


C{u/v) — filogn 
oVlogii 


<t = 


\/27T 7-0 


2 dx, 


and similarly for P„ and Tl n with the same constants. The mean fi satisfies 


(1) 


H := 


m 


E c ( n ) / , €{x)dx 


where f )(T) := log | T' (x)\f(x)dx — 6 f og2 - If we define if(x) := c(n) + fi log x for 

all x £ (—^r, i] and i[>(0) := 0, then the variance o 2 satisfies 


a 2 = 


m 


lim — 

n—>oo Tl 


j o ^E^*^ > 0 . 


As an example of the application of Theorem [T] if c is set to be the constant 
function 1 then the cost C(u/v ) corresponds to the number of division steps required 
to evaluate gcd(u,n). Alternatively, if we define c(m) := 1 and c(n) := 0 otherwise, 
then the resulting cost C(u/v) is the number of occurences of the digit m in the 
continued fraction representation of the number u/v. 

By a classical theorem of V. A. Rokhlin m the quantity \)(T) is precisely the 
entropy of the transformation T with respect to the absolutely continuous invariant 
probability measure on [0,1] with density However, this identification has no 
bearing on our argument beyond its influence on our choice of notation. While our 
exposition in this article restricts itself to the classical Euclidean algorithm, it is 
easy to modify our method along the lines of [TJ .25] so as to apply to its Odd and 
Centred variants. 

The proof of Theorem |T] which we present uses the method of moments. By 
standard results of probability theory, to prove that the distribution of C(u/v) 
converges to a normal distribution it is sufficient to prove that for every integer 
p > 0 the sequence of centred p th moments, 


( 2 ) 


1 


E 

(u,v)£Q,n 


C(u/v) — log n 
oVlogn 


converges to the p th moment of a standard normal distribution (see e.g. [2} §30]). 
To study these moments we consider the bivariate Dirichlet series 

oo 1 

:= E“2l E exp(u(C(u/v)~ p,\ogv)). 

n=l 


The sum which appears in the p th moment © resembles the sum of the first n 
terms of the p th partial derivative with respect to ui of D(s,lo) evaluated at s = 
lo = 0, differing in the presence of filogn versus fi log v inside the summation. 
To exploit this fact we extract this sum from D(s,w) using a Tauberian theorem 
due to H. Delange, stated below as Theorem [2] When p is even Delange’s result 
may be applied directly to the calculation of the moment. For odd p our analysis is 
complicated by the fact that the summands fail to be non-negative, and to overcome 
this issue we adopt a strategy due to H.-K. Hwang and S. Janson mm - These 
parts of the proof constitute ]J5] below. 
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To show that the Dirichlet series has the properties needed to apply Delange’s 
theorem we equate the series with a summation over all of the possible compositions 
of inverse branches of T. This is achieved by expressing these sums over inverse 
branches in terms of the Gauss-Kuzmin-Wirsing operator, the necessary features 
of which are studied in [J3] below. In 2] we effect the translation of the properties 
of the transfer operator into those of the Dirichlet series. 

In preparing this article we have taken some effort to minimise the prior knowl¬ 
edge of ergodic theory and thermodynamic formalism required on the part of the 
reader. In particular, the only significant “black box” results to which we appeal are 
Delange’s Tauberian theorem, the efficacy of the method of moments, and certain 
results from the perturbation theory of compact operators, for the last of which we 
refer the reader to the book by T. Kato [Hj. The decision to restrict our attention 
to bounded costs is partly based in this effort at conciseness: to treat the general 
case would require either appeal to an external result in thermodynamic formalism 
(as occurs in CP) or an exposition which would have substantially lengthened (}3l 


3. Part I: the transfer operator 


In this section we define a weighted version of the Gauss-Kuzmin-Wirsing oper¬ 
ator and establish those spectral properties of that operator which will be used in 
the subsequent sections. This material is to a significant extent expository and is 
included for the sake of completeness. 

Here and throughout the article we let D:= {zgC: |z— 1 | < 1 }) and let £: D — > 
C be the holomorphic extension of the function £ defined in the introduction. Since 
the function z4z + n has no zeros in ED we may unambiguously write (z + n) s = 
exp(s log(z + n)) for every zeD and s £ C, where log(z + n) indicates the branch 
of the complex logarithm which coincides with the real logarithm when z £ [0,1]. 
We recall that a function from an open subset of C to a Banach space X is called 
holomorphic if it is locally equal to a convergent power series with coefficients 
in X. A function from an open subset of C to X is holomorphic if and only if 
its composition with each element of X* is holomorphic in the ordinary sense. 
Throughout this article we shall repeatedly appeal without comment to the fact 
that a uniform limit of holomorphic functions is holomorphic, which follows from 
the combination of Cauchy and Morera’s theorems. We recall that H°°(] D>) denotes 
the set of all bounded holomorphic functions on B, which by the aforementioned 
principle is a complex Banach space when equipped with the uniform norm. 

Here and throughout the article we let |c|oo denote the least upper bound of the 
cost function c, and let fj, > 0 be as defined by ©■ For each s, u £ C with 5R(s) > 4, 
and each / £ ff°°(B), let us define a pair of functions £ S)W /, ^s,wf '■ B —> C by 



(3) 



(The difference between these two functions - that the sum defining J- s ^ does not 
include the index n = 1 - reflects the fact that the continued fraction expansion of 
a rational number is defined unambiguously only if the final digit is not permitted 


to equal 1.) We begin the proof of Theorem [T] by quickly proving the following 
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result which is of a somewhat classic type; compare for example [16] Prop.l], [25j 
Thm. 1], ;T| Prop. 0]. 

Proposition 3.1. The equation J3]) defines a family of linear operators C s T StU1 
on H°°(D) such that: 

(i) For every {s,lo) £ C with 5ft(s) > ^ and |w| sufficiently small, C S}U1 and 
Fs,ui ore well-defined compact operators on H°°( D). The dependence of these 
operators on (s, w) is holomorphic. 

(ii) The operator Cip acts boundedly on L 1 ([0,1]) and satisfies fg (£i,of)(x)dx = 
fg f(x)dx for all f £ L 1 ([0,1]). If f,g: [0,1] —> C are such that f and 
f(goT) both belong to L 1 ([0,1]), then Ci : o(f(g°T)) = (Ci t of)g Lebesgue 
almost everywhere on [0,1]. If f € i 1 ([0,1]) then also f o T £ L 1 ([0,1]) and 
\f °T\l^ < C|/| L i where C > 1 is constant. 

(Hi) There exists an open set V C C 2 containing (1,0) such that for all (s,w) £ V 
we may write C SiUJ = X(s,vj)V s ,u + Af s ,u where V s ,uA f s , u = A f s , u V s ,uj = 0, 
V s is a projection onto a one-dimensional subspace of H°°(D) which is an 
eigenspace of C SiUJ , M s , u has spectral radius strictly less than 1, and X is a 
complex number. The dependence of X,V and J\f on (s,w) is holomorphic 
throughout V. We have A(1,0) = 1 and V\pf = (fg f(x)dxj ( for all f £ 
fP°(B). j 


Proof, (i). For each n > 1 and / £ H°°( B) it is easily seen that the equation 


(B a , u , n f) (z) := 


0 uc{n) 


if 


(n + z) 2s+flul \n + Z/ 

defines a holomorphic function on B. Since \z u \ = e~ ar s( z )when z,u £ C 
with > 0, we furthermore have 

2 


(4) 


\l3 s ,uj,nf\oo < 


Jc|ooSf(w)+7r(a(s)+f Ss(ui)) | 


In particular B s ,ui,n is a bounded linear operator on 7J°°(B), and it is clear that 
Bs,ui,n depends holomorphically on (s,w). Since C S)U] is simply the sum over all 
n > 1 of the operators B s ^^ n , and for JJ(s) > 4 and |w| small enough the inequality 
0 implies that the convergence of this sum is locally uniform with respect to (s, w), 
C s ^ u is a holomorphic family of bounded operators on .ff°°(B). Taking the sum over 
n > 2 shows the family F StU] to have the same properties. 

To complete the proof of (i) it suffices to show that each operator B s ,ui,n is 
compact. If h: B — > B is any of the linear fractional transformations z i —> 1 /(z + n) 
then we may expand its domain of definition to the larger disc B := {z £ C: \z— 1| < 
4^}, and note that h(B) is contained in a compact subset of B. The elements of the 
set {foh: f £ H°°( B), |/|oo < 1} may therefore be simultaneously extended to the 
larger domain B whilst remaining uniformly bounded and holomorphic. It follows 
by Montel’s compactness principle that this set is a precompact subset of H°°( B) 
and hence the linear operator / i —> f o h on H°°( B) is compact. If g £ H°°( B) 
is any bounded function then the map / i —> g ■ (/ o h) is the composition of the 
compact operator / i —> f o h with the bounded operator f g ■ f and thus is a 
compact operator on U°°(B) with norm |g|oo- By this principle it follows that each 
of the operators B s ^,n is compact. 
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(ii). For all / £ i 1 ([0,1]), by taking the substitution y = 1 /(n + x) we obtain 


<£,,„/) (x)dx = l e J dx = £ 79 = l 1{x)± 


n= 1 n + 1 

and an obvious modification shows that \Ciflf\L 1 < l/U 1 - If / £ L 1 ([0,1]), 
g: [0,1] —> C is measurable and x £ [0,1] then 


(£i,o(/-(9or)))W = E(iV3/ 

n=l ' ' 


l 


x + n 


9 T 


1 


= ffM E 


f 


—' (x + n ) 2 ' \x + n 


n + x y 

= 9{x)(Ci,of)(x) 


so that Ci : o((g ° T)/) = g(C\fif) almost everywhere as claimed. In particular if 
/ £ F 1 ([0,1]) is essentially bounded then clearly / o T £ L 1 ([0,1]) and we obtain 
fo \f(Tx)\dx = fg £i, 0 (|/ ° T\)(x)dx < |£i, 0 l|oo fo \f(x)\dx, which proves the last 
part of (ii) for /; the general case follows by approximation. 

(iii). Since £ 1 ^ is compact the nonzero elements of its spectrum are all isolated 
eigenvalues of finite multiplicity. By direct calculation £i,o£ = £■ Let £ be an eigen¬ 
function of £ 1,0 corresponding to some eigenvalue A, and by rescaling if necessary 
suppose that sup xe r 0 ^ |£(x)|/£(x) = 1 with this supremum attained at xq, say. We 
have 


|A^(ic 0 )| = |£i,ol(zo)| = 


y - £ 

“ {xq +n) 2 \xo + n 

n= 1 


^E 

n= 1 


(x 0 + n) 2 


1 


xq +n 


^E 


(xq + n ) 2 \xq + n 


= (£ 1 , 0 ?) (x 0 ) = £(x 0 ) = |l(zo)| 


and therefore |A| < 1, so £ 1,0 has spectral radius 1. If |A| = 1 then the above 
inequalities must be equations. Since |£(x)| < £(x) for all x £ [0,1] we deduce 
that |£(l/(xo + ra))| = £(l/(xo + n)) for all n > 1, since otherwise the second 
inequality would be strict; and moreover the argument of £(l/(xo + n)) must be 
independent of n since otherwise the first inequality would be strict. It follows that 
£(l/(xo + n))£(l/(xo + ?r)) _1 is constant in n, and since l/(xo + n) —> 0 £ D as 
n —> 00 we conclude that the holomorphic function £/£ must be constant on D. 
We have thus seen that ker(£io — Id#oo( D )) is one-dimensional, that p(£i,o) = 1, 
and that £±.o has no other eigenvalues of modulus one. To see that £ 1.0 has no 
generalised eigenfunctions at 1 we observe that if £ £ ker(£i >0 — Id#oo( D )) 2 then 
(£ 1,0 — Idff°°(B))£ must be proportional to £, but by part (ii) this function has zero 
integral on [0,1] and hence must be zero. We conclude that ker(£i i o — Id#oo( D )) 2 = 
ker(£i ; o — Id#o°( D )) and 1 is a simple eigenvalue of £ 1,0 acting on H°°( D). 

Let r C C be an anticlockwise-oriented closed curve which encloses 1 but does 
not enclose any other point of the spectrum of £ 1 , 0 - By m p.212] we may find an 
open ball V containing (1, 0) such that the spectrum of £ SjW does not intersect T for 
any (s, w) £ V. Let V s ,u> := — ^ / r (.z-Id#oo( D ) — £ SjW )" 1 dz be the R.iesz projection 
associated to the curve T, which is a projection commuting with £ S)W on p-178]) 
and depends holomorphically on (s,u>) (in p.369]). The image of V 1.0 is a closed 
subspace since that operator is a continuous projection, and is preserved by £po 
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since that operator commutes with V\,g. The spectrum of £ 1,0 restricted to the 
image of Vi,o is precisely {1} (by [12 P-178]) and in particular the restriction of £yo 
to that space is both compact and invertible. The image of Vi,g is therefore finite¬ 
dimensional, and since 1 is a simple eigenvalue of £ 1.0 it follows that rank^yo = 1- 
By P2J p.212] we thus have rank7^,^ = 1 for all (s,w) G V. The function V s ,uf, is 
thus an eigenfunction of £ s , w corresponding to some eigenvalue A(s, w) G C. Define 
A F s< w := £ s ,^(Id ff c=o( D ) - V a< u) to obtain C s>u = X(s,w)P a>u + Af s , u - The spectrum 
of A/yo equals the spectrum of £yo with 1 removed (pdj p.178]) so the spectral 
radius of A/ 1,0 is less than one, and by upper semicontinuity of the spectral radius 
we may replace V with a smaller set to obtain p(Af s , u ) < 1 for all (s,co) G V. 

For each / G H°°(l 0) we have jCi,gPi,gf = 'Pi,of so that Pi,of is proportional 
to £, and since p(A/ 1 , 0 ) < 1 we have f* (Pi,of)(x)dx = liinn^oo fg(jC™ 0 f)(x)dx = 
Jo f( x )d x using (ii), which shows that P\,gf = (fg f(x)dx)t ; as required. Since in 
particular Pi ,oC = £ it follows that (Ps,u>€)(0) / 0 when (s, to) is close to zero, so by 
taking V smaller if necessary we may write A(s,u>) = ((P s ,u J f)(0))~ 1 (Ts,ujPs,ujO(0) 
for all (s,co) G V which shows that A is holomorphic. □ 

The elimination of unit eigenvalues on the line 5ft(s) = 1 is classical when transfer 
operators interact with Tauberian arguments, see for example m Prop. 6.1]. For 
versions of this argument in the context of transfer operators associated to number- 
theoretic algorithms, see for example [13 Lemma 6.15],[25] §5.5]. 

Proposition 3.2. Let s G C such that 5ft(s) > 1 and s ^ 1. Then p{C s . 0 ) < 1. 

Proof. Let us first show that if 5ft(s) > 1 then p(C s ,o) < 1 - Using Proposition 
Eli) we may choose an eigenfunction of C s ,g which corresponds an eigen¬ 
value of modulus p(C s , 0 ). By rescaling this function if necessary we assume that 
sup xg r 0 i] |£ s (a;)|/£(x) = 1 with this supremum attained at some point xo G [0,1]. 
For every n > 1 we have 0 < (xg+n)~ 2 < 1 and therefore |(a;o+n) _2s | < {xg+n)~ 2 . 


Thus 




and therefore p(C s ,g) < 1 as claimed. 

We now vary this argument to show that p(Ci+u,g) < 1 for all nonzero real 
numbers t. Given t G R, let be an eigenfunction of £i+it,o which corresponds 
to an eigenvalue of modulus p(Ci+it,g). In the same manner we rescale so that 
sup^grQ |£t(x)|/£(2;) = 1 and choose xg G [0,1] attaining this supremum. We 
similarly estimate 



(£i,oO ( x o) = f.{xo) = |6(a>o)| 
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which in particular establishes p(C\+it,a) < 1. Let us suppose that p{Ci+it,o) = 1, in 
which case both inequalities above are equations. As with the proof of Proposition 
cnii) this is only possible if the argument of (xo + n) 2 2,t ^t( l/(xo + ra)) is inde¬ 
pendent of n and the modulus of that expression is equal to (xo + n) _ 2 £(l/(xo-|-n)). 
It follows that for some fixed 6 G R and for all n> 1 


je 


(xq + n) 2+2it ' \xq + 


(x 0 + n ) 2 \ Xq “t“ 


and therefore 


lim 

n—>-oo 


e 


iO 


(xq + n) 2lt 



m 


but the former sequence is divergent if t 7 ^ 0 . We conclude that if p(Ci +itj o) 
then necessarily t = 0 , and this completes the proof of the proposition. 


= 1 
□ 


Proposition 3.3. Let ip: [0,1] —> C be as defined in Theorem 
satisfies 


— (l, 0 ) = -f,(T)< 0 , 


53 (1 '°) = 0 - 


0 


The function A 





Proof. Let us write := V s ,w£, when (s,u>) G V so that = A(s,w)£ SitJ in 

this region. Since Vigt; = £ we have £ 1,0 = £. Since V s)W is holomorphic, is 
infinitely differentiable in H°°(W) and hence also in L 1 ([0,1]): we use the notation 
to refer to the first partial derivatives with respect to s and u, and for 
the second derivative with respect to ui. We denote the corresponding derivatives 
of A by A S ,A W and A ww . 

We claim that s 1 —> |T , | 1_s £ s ,o is a holomorphic map from a small neighbourhood 
of 1 into L 1 ([0,1]). By the holomorphicity of s 1 —> £ Si o near s = 1 we may write 
£s,o = S^=o( s ~ l) n 5n where each g n G and limsup^^ \g n \\L n < 00 . On 

the other hand |T'(x)| 1_s = x 2_2s = x ) n (2 ~ 2s) n /n! for x G (0,1], and 

since fg | logx|"dx = / 0 °° y n e~ v dy = n\ it follows that the function s i-A |T'| 1_S is 
a holomorphic map from a neighbourhood of 1 into L 1 ([0,1]). It is now a simple 
exercise to write s 1 —> |L 1, | 1_s £ s ,o as a power series in (s — 1 ) with coefficients in 
L 1 ([0,1]) and a nonzero radius of convergence, which proves the claim. For s close 
to 1 and x G (0,1) we may write 


(£s,o6,o)(z) = 


1 


= £ 


(x + n) 

1 


(x + n) 


2s ^ s ’ ( 

't' 


1 


x + n 

1 

x + n 


£s,0 


1 


x + n 


= A,o (ITf-^sfi) (x) 


and hence the equation C\g (|X 1 '| 1 s £ s ,o) = A(s,0)£ Sj o is valid in L 1 ([0, 1 ]). Differ¬ 
entiating this equation at s = 1 we obtain 

A,0 (- log |T'|£ + £,„) = A s (l, 0 )£ + £1% 

since £ 1,0 = £. Integration yields A S (1,0) = — fg log |T'(x)|£(x)dx = — f)(T) via 
Proposition 13. lf ii) as desired. 
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We next claim that for each m > 1 the map to H > exp(w V’ ° T k )£i tU is a 

holomorphic function from a small neighbourhood of 0 G C into L 1 ([0,1]). Since 
we may write for some sequence of functions h n G 

with limsup^^go |/i n |oo < oo, it suffices to show that to K > exp(w Y^k=o ^ ° T k ) 
is a holomorphic map into L 1 ([0,1]). Since |^>(x)| < |c|oo + |2^ilogx| for all 
x G (0,1] it is clear that J 1 \ijo{x)\ n dx < M n n\ for some constant M > 0 de¬ 
pending on c, and by Proposition IQ.ll' iil there is a constant K > 0 such that 
fg \il)(T k x)\ n dx < K k fg \ij){x)\ n dx < K k M n n\ for every k > 0. Using the gener¬ 
alised Holder inequality 


m— 1 


E ^ Tkx ) 


k—0 


m— 1 i n— 1 

dx < e / n dx 

fci,...,fc„= 0 "'° i =0 

m—1 n—1 / 

< E n / \^{T ki x)\ n dx) <m n M n K m n\ 

fci,...,fc„=0 i=0 w° 


and it follows that the power series ^ (Efclo 1 V’ ° / n - with coefficients 

in Zr([0,1]) has nonzero radius of convergence, completing the proof of the claim. 

Now, if to is close to zero then since A(l,w)£i iaj = £i, w £i,uj, using the definitions 
of C\ tUJ and ij) we have 


(5) A(l,w)£i iW (x) = E 7T 

n =1 v 


+ nf 




= £ 


1,0 


UJ'lfj 


6,w) (x) 


for almost every x G [0,1]. We claim that for all n > 0 the equation 

(6) £i, 0 (e“ *° Tk &,„) = A(1, U )% 

in U 1 ([0,1]) is satisfied for all to sufficiently close to zero. Clearly the claim holds 
for n = 0. Given the validity of the claim for some n > 0, if to is small enough that 
additionally e“^£=o ^° Tk £i,u> € i 1 ([0,1]) then we may deduce 

£«+! ( e "£2=o^) =Cl 0 {e u ^=o^ oTk C lfi (e^^)) 

= A(l,c o)C n lfi = X{ i^)n+i^ 

using Proposition 13.11 ii), (O and ©, and the claim follows by induction. Differen¬ 
tiating © in U 1 ([0,1]) yields 


£"o U E *=° V ' oTfc (EE ° T k )^ + Sk 

V Vfc=o 

= nA(l,cur- 1 A w (l, W )ei, w + A(l, W )”^- 

Setting n = 1, to = 0 and integrating over [0,1] gives A u (l, 0) = fg ip(x)f;(x)dx = 0 
as claimed. Differentiating a second time, setting to = 0 and substituting £ 1,0 = £ 
and A w (l, 0) = 0 yields the expression 



£ 


n 

1,0 





= nA* w (l,0)£ + £$. 
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Integrating, cancelling / 0 ' (x)dx from both sides and dividing by n we find that 

A™(1,0) = i jf' £(*) + 2 dx 

for all n > 1. Using Proposition 13.lf iiil we have 

i n —l n —1 / -l \ 

linr -^^1,0= -E^“o= / ■ € 

n—»oo Tl z —' ’ n—>• oo 77 , z —' ’ ’ \ /n / 

fc=0 /c—0 WU 7 

in the uniform norm, and therefore 

-I pi n -1 pi / n—1 \ 

(7) linr-/ ^''ip(T k x)^i 0 (x)dx= lim / i/-^) - ^ 4,o^,o W* = 0 
n ^°° n J o ^ ^°°7o V 71 ^ / 

using Proposition 13.lf iii on each summand together with f n if>(x)t;(x)dx = 0. Thus 
A WW (1)0)= lim - [ fy^ i/j(T k x)] £,(x)dx> 0 

"^°° n ^ \k=0 J 

as required. It remains to prove that this quantity is nonzero. Define 

^ c(n) - filogjz + n) ^ ( 1 


(z + n) 2 


z + n 


so that T G H°°( D), and note that £i,o(^) = ^ almost everywhere in [0,1], so 
in particular ty(x)dx = fg ip(x)£(x)dx = 0. It follows that 
for every n > 1 and hence the sum x ’■= £ _1 (£i o^O defines an element of 

H°°(V>). Define J) := 'ip + x ~ X° T G T 1 ([0,1]) and observe that = 

dt + £i,o(x£) — = 0 in T 1 ([0,1]). For each n > 1 we have 

/ n—1 \ 2 / n—1 \ 2 /n—1 \ 

E ^ ° - E^ oTfc =2(x°T n -x) E^ oTfc +(xoT”-x) 2 


\k =0 


\k—0 


\k—0 


in L 1 . By an argument similar to 0 we have A JJ'k=o 4’°T k —> 0 in T 1 ([0,1]), and 
since x ° T N — y is bounded on [0,1] independently of n we may easily derive 

2 


( 8 ) 


A WW (1,0) = lim 

n—too 


\j Q ^E^( T E)j Z(x)dx. 


Now, using Proposition 13.lf iii and the identity £i,o£ = £, for k > j > 0 we have 
J f ip(T k x)ip(T j x)Z(x)dx = j (j) (T k x) ip(T j x)t(x)^ dx 

= [ ip(x)C k ~ 3 (ipO(x)dx = 0 
Jo 

since £i,o(V’0 = 0. The off-diagonal terms of the squared sum in 0 thus vanish, 
leaving 

Awo)( 1 5 0)= lim - [ ^(T k x) 2 j £(x)dx = f ■ip(x) 2 t;(x)dx. 

n->oo n J Q J J 0 
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If this quantity were zero then we would have ip = %oT —x almost everywhere by the 
definition of ip. Since c is bounded and the logarithm function is not, the function 
ip is not essentially bounded; but x is bounded, so the equation ip = % o T — a.e. 
is impossible and we conclude that A WW (1,0) > 0 as claimed. □ 


4. Part II: the Dirichlet series 


This short section comprises the following result: 


Proposition 4.1. For each integer p > 0, the Dirichlet series 

oo 1 

D p( S ) : =^2~ 2 l E (C(u/v)-p \ogv) p 

n— 1 (u,i;)£Q n \Q n _ i 


is absolutely convergent in the half-plane 5ft(s) > 1. For each odd p, D p (s) is 
holomorphic on the line 5ft(s) = 1 with the possible exception of a pole at s = 1 of 
order not greater than (p+ 1)/2. For each even p there exist a holomorphic function 
R defined in the half-plane Sft(s) > 1 and complex numbers «i,..., a p / 2 such that 


D P {s) 


P ] -cr p fo 

2Pt){T){s- l) 1+ 5 


£ 

2 

+E 

k= 1 


dk 

(s - l) fe 


R(s) 


when 5ft(s) > 1, where fo := (T- i,o£i,ol)(0) > 0 and o 2 := 2X bJbJ (l, 0)/[)(T) > 0. 


Proof. For each n > 0 let us define S n to be the set of all u/v £ (0,1) which have a 
continued fraction expansion u/v = [a\,. .. , a n + 1 ] with a n + 1 > 2; equivalently, if we 
write h n (x) := 1 /{x+n) then u/r e S„ if and only \iu/v = (h ai oh a2 o- ■ -oh an+1 )(0) 
with a n +1 > 2. This composition of functions hk is itself a function of the form 
x i ^ (ax + (3)/(xx + 5) where \a8 — /3q| = 1 and a, (3,^,5 £ Z, so its derivative 
at x = 0 is equal to 1 jo 1 = 1/v 2 . Since each u/v has a unique continued fraction 
expansion such that a n +i > 2, using the definition of F S)UJ and C s>u] we may write 

(J-,,^1) (0) = Y, 

a\ ,...,a n >l 

CLn -\-1 >2 

= E ^2iexp(w(C'(w/t;)-/ilogi;)) 

(u,v)ES n 

for each n > 0. Notice that when oj and s are real this is a series of non-negative 
real numbers whose total is bounded by |.Fs,w|oo|£g w |oo- Taking the sum over all 
integers n > 0 we obtain 

OO -J OO 

E^ 2 l E exp (w (C(u/u) - plogv)) = Y (Fs.uCs^ 1 ) (°)- 

n=i fc=o 

The second of these infinite series is well-defined for all (s, v) such that p(C StUJ ) < 1. 
The first series is a convergent series of non-negative real numbers when (s, w) are 
real and p(£ s ,w) < 1, and it is not difficult to see that this implies the absolute 
convergence of the series when p (£sr(s),!R(lj)) < 1- By Proposition 13.21 both series 
converge absolutely and m is a valid equation when 5ft(s) > 1 and w is sufficiently 


(h ai o ■ ■ ■ o h an+1 )' (0) 
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close to zero. When these conditions are met we may differentiate both sides with 
respect to ui to obtain for each integer p > 0 
( 10 ) 


n =1 (ii,u)Gr2 n \f2 n _i 


g P 

dujp 


OO 

E (^(.•)4,-) 1 ) (o) 


M) 


with the former series being absolutely convergent when Jf(s) >1. In view of 
Proposition l3.2l the second series in © also convergent when 3?(s) = 1, s ^ 1 and w 
is sufficiently small, and so the right-hand side of (fTOl) is holomorphic in the region 
9?(s) > 1, s / 1. In particular the series D p (s ) admits an analytic continuation into 
that region as required. 

The remainder of the proof is concerned with the analysis of D p (s ) close to s = 1. 
Let us fix the integer p > 0. Using Proposition 13.11 iiil we may find an open set 
V C C 2 containing (1,0) in which C s<u] = A (s,uj)V StUJ ©A f a>u where A(1,0) = 1, 
V s ,ui is a projection of rank one and the spectral radius of A/), w is strictly less than 
1. When (s,w) £ V, !fi(s) > 1 and uj is sufficiently small we have |A(s,w)| < 1 by 
Proposition 13.21 and in this case 

OO OO 

E (0) = E A(s,^) fc (Fs^V^l) (0) + (0) 

0 k=0 

= 1- 77 - +N(s,U>), 

1 — A (s, uj) 

say, where N is holomorphic throughout V. For (s,0) £ V with 5f(s) > 1 the 
left-hand side of C© is therefore equal to 

(0)) M ) + Ri(s) 

where R\ extends holomorphically to a neighbourhood of s = 1. Let us expand the 
first expression as 


E 


n=0 


n 

p 


(rh) ( *’°>) m) - 


Since (J r SjW 'P SiW l)(0) is holomorphic throughout V its partial derivatives are holo¬ 
morphic also, and hence poles may only enter the above expression via the partial 
derivatives of 1/(1 — A). We will show that the term corresponding to n = p con¬ 
tributes a pole of order [1 +p/ 2J, and the remaining terms can only contribute 
poles of lower order. Using Faa di Bruno’s formula (see e.g. [20]). 

( 11 ) 

d n ( 1 \ ^)- tt AL^(5,0) fc< 

^ Vi - aE ’ j fcl+2fe2 EWn (1 - A(-,0))^*« 11 km k < 

k lt ...,k n >0 


We consider separately the contribution from each summand. 

Recall from Proposition l3.3l that A w (s, 0) has a zero at s = 1 with degree equal to 
some positive integer d , whereas 1 — A(S, 0) has a simple zero at s = 1. If k\ > 2 then 
X )"=2 ^ — \( n — ^i) £ n /2 — 1 and so the function A w (s, 0) fcl (l — A(s, 0)) _1 “^ ke 
has a pole of order at most n/2 at s = 1. If ki > 1 for some i > 3, then X )"=2 ^ — 
1 (n — hi) < (n — l)/2 and hence A w (s, 0) fcl (l — A(s, O))” 1 ”^ k e has a pole of order 
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not greater than (n+ 1)/2. The only remaining case is that in which k 2 = \n/2\ and 
k\ = n— 2 k 2 , and in this case the order of the pole of \ u (s, 0) fcl (1 —A(s, 0)) -1- ^ kl 
at s = 1 is precisely 1+ [n/2j +(1—d)/ci. It follows that when n is odd the expression 
CD has a pole at s = 1 of order not greater than (n + l)/2, and when n is even it 
has a pole of order precisely 1 + n /2 at s = 1 which arises solely from the summand 
with k 2 = n/2 and k\ = 0. Summing dill) over n in the range 0 to p , it follows that 
whenp is odd D p has a pole of order at most (p+l)/2ats = l, which completes the 
proof of the proposition for odd p. In the case where p is even the above analysis 
shows that we may write 


d p 

dujP 


1-A 


OF (v) P ( . >0 l)(0) M) 


p!A^(s,0) p/2 
2p/ 2 (1 - A(s, 0)) 1+p / 2 


s, 0 v a , 0 im + R2 (s) 


for all s / 1 such that (s,0) G V, where R 2 has a pole at s = 1 of order not 
greater than p/2 and otherwise is holomorphic at s for ( s , 0) € V. Let fo := 
(J r i ! o7 : ’i ! ol)(0), which is positive since V\pl = £ and preserves the set of 
functions which are positive on [0,1]. Since A(1,0) = 1, A WW (1,0) = \ c 2 h(T), and 
|A(s,0)| < 1 when SR(s) > 1 and s / 1, the functions ^4(s) := (1 — s)/(l — A(s,0)), 
B(a) := (A ww (s, 0) — ^a' 2 t)(T ))/(1 — s) and C(s) := ((^,o^ s ,o(l)(0)-fo)/(l-s) are 
holomorphic in the region 9?(s) > 1. Substituting (1 — A(s,0)) _1 = (1 — s) -1 A(s), 
A ww (s, 0) = icr 2 t )(T) + (1 - s)B(s) and (J' s ,o'P 5 ,ol)(0) = fo + (1 - s)C(s) into the 
previous expression we obtain 


£ M) = 


p\aP\){T)P/ 2 A{s) 1+p / 2 fo 


R 3 (s) 


2p(1 - s) 1+ p/ 2 

for all s sufficiently close to 1, where R 3 extends holomorphically to the region 
5R(s) > 1 with the possible exception of a pole at s = 1 with order not greater than 
p/2. In view of Proposition EHfl we have Al(l) = A S (1,0) -1 = — l/b(T), and so this 
may be further rewritten as 


d p 


1 


)M) = 


p!er p fo 


2Pfj(r)(s- iy+p / 2 


I?4(i 


where R 4 has the same properties as R 3 . The proof is complete. 


□ 


5. Part III: the Tauberian argument 


In this section we apply the properties of the Dirichlet series D p (s ) established 
in Proposition 14. 1 1 to prove Theorem [TJ We give the full details of the proof only in 
the case of the distribution of costs on O n ; the case of the distribution on fl n can 
be handled by simple modifications which are indicated at the end of the section. 
By the classical method of moments, to prove Theorem [1] it is sufficient to prove 
that for every p > 1 the sequence of p th moments 


converges to the limit 


1 


E 


(u,v)€Q„ 


1 



t p e 2 ^ dt = 


C(u/v) — plogn 
\/log n 


p 


0 

<r p (p — 1 )!! 


if p is odd 
if p is even, 


see for example j2] §30]. We will begin by calculating the asymptotics of the related 
sequence ^ 2 ( u (C( u / V ) ~ d^ogv) p as n —> 00 for each p > 0, and then use this 
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to derive the corresponding asymptotics for the sequences of moments. We require 
the following Tauberian theorem due to H. Delange [6], Th. Ill]: 


Theorem 2 (Delange). Let (a n ) n >\ be a sequence of non-negative real numbers. 
Suppose that for all s in the half-plane 5ft(s) > a the series D(s) := a nn~ s is 

absolutely convergent and satisfies 


D(s) 


g(«) t y- 9j(s) 

(s — a) u 2 —'(s — a) x j 


+ h{s) 


where w > Xj > 0 for every j, the functions gi, g and h are holomorphic in the 
region 5ft(s) > a, and g(a) > 0. Then 


( 12 ) 


lim _I_ y a - g(a) . 

n^-oc N a (\ogN) ul 1 " ar(w) 


Remark. Delange’s statement includes some additional parameters fij all of 
which we set to 0 in the above statement. Since we will apply the result to series 
of the form a nn~ 2s we obtain an additional factor 2“ _1 in the numerator of 

the limit in (1121) which arises from the change of variable. 

Let us write D^ for the fc th derivative of D p with respect to s. In the case 
p = 0 the combination of Proposition lTTI with Theorem [2] yields linin^oo ffLl n /n 2 = 
fo/fl(T). For even p > 0 the same approach yields 


lim — 

n —>oo 77, z 


1 

(logn) p / 2 


E 


(C(u/v) 


p log v) p 


P } -v p fo 

2P/ 2 (p/2)!h(T) 


and therefore 

(13) lim ---- —t- 

n^oo fffl n (\ogn ) p / 2 


E {C{u/v)~ g\ogvf = 


p\a p 


(u,v)(£Q n 


2P/ 2 (p/2)\ 


= (p — l)!!cr p . 


When p > 1 is odd the coefficients of the series D p fail to be non-negative and a 
more delicate argument is required, which comprises much of the remainder of this 
section. We adopt an argument due to H.-K. Hwang and S. Janson mm. Let us 
fix odd p > 1 and define auxiliary Dirichlet series, 

- Eji £ ((c0 - flog.) 2p +(log.)”). 

n=1 (u,v)eQ n \Qn-i ^ ' 

00 1 2 
V ^ s ) : = E^I E - M logu) P + (logi;)5) , 

n= 1 (u,v)£Qn\Qn— 1 

OO 1 

(14) V 3 (s) ■= E (!ogu) f (C0 - m log P , 

n= 1 (u : v)£Q n \Q n _i 

where we note that V 3 = \{T> 2 — T>i). The first of these series is easily seen to satisfy 
T>i (s) = D 2p (s) — j^Do P \s), and since D 2p and D 0 converge absolutely in the region 
5ft(s) > 1 the series T >i converges absolutely in that region also. Since the terms 
of V 3 (s) are bounded in absolute value by those of (— 2) _ ^p/ 2 1 (Dp^ p / 2 ^(5i(s)) + 
Dq P ^ 2 1 (JJ(s))) it follows that V 3 (s) converges absolutely for 5i(s) > 1, and we 
deduce that V 2 = V\ + 2V 3 has the same property. 
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Using the properties of D 2 P and Do described in Proposition 14 .1 1 we may write 

g{s) 


■DM = D 2r (,) - ±D<r>(.) = (t _ 1)I+p 


/(s) 


where g and / are holomorphic in the region ?R(s) > 1 and g( 1 ) is positive (its 
precise value is unimportant). Using Theorem [2] it follows that 

2 p 3(l) 


lim 


1 


n -¥oo n 2 (log n) p 


^2 ( (C(u/v ) - \i log v) 2p + (log v) p ^j 




Pi 


We will show that T >2 has a singularity at s = 1 with the same order and leading 
term as T>\ and hence that its partial sums have identical asymptotic behaviour. 
To achieve this we shall bound the singularity of the difference T> 3 . We have 



D 


m 


» = £ 


n= 1 


- y 

i2s / j 

(u,v)&Q n \Q n _ 1 


(log»)^(c0 


- M 



P 


for 5ft(s) > 1. Noting that v = n when (u,v) £ Cl n \ U ra _i and using the identity 
/ 0 °° n~ 2 t t~^dt = i/ 7 t /(2 logra) we may rewrite (fill) as 

OO 1 

y^ y (logc) 2 ^ (c(^)-/iio g u) p 

n =1 

(-l) 2 ^ f°° D { p {p+1)/2) (s + t) ^ 

V2^ Jo Vt 

We have \Dp P+ 1 ^ 2 \s)\ = 0( 2 -K ( 2s )) uniformly in S(s) in the limit as 5ft(s) —> 00 
since Dp P+1 ^ 2 ' > is a Diriclilet series with zero constant term. It follows that when 
Dp( p+1 )/ 2 ) is bounded and holomorphic in a strip of the form {s: 5ft(s) > a and /? < 
3(s) < 7 } the integral / 0 °° Dp P+1 ^ 2 \s + converges absolutely uniformly 

in the same strip. We deduce that V 3 is holomorphic in the region Ji(s) > 1 except 
for the point s = 1 . 

By Proposition 14.11 l) r has a pole at s = 1 with order at most [p + l)/2, and 
therefore D$ p+1 ^ 2 '* has a pole at s = 1 of order at most p + 1. It follows that we 
may find an open set U of the form U = {s: Sft(s) > 1 — e, |3(s)| < 5}, a bounded 
holomorphic function R defined on U , and complex numbers a±,..., a,( such that 


dt 


0 n 


~ 2t Vt 



Dp (p+1 )/ 2 ) (s) = y 


CLk 


k= 1 




R(s) 


for all s £ U, where t < p+ 1. We have |i?(s)| = 0(5ft(s) 1 ) as Sft(s) 00 uniformly 
throughout U in view of the corresponding bound on |D^ p+ 1 ^ 2 ^(s)| and hence the 
integral J 0 °° R(s + f)f _1 / 2 df converges absolutely uniformly in this region. Since 

r 1 -k r du n(2k — 2)! ±_ h 

1 - 7 dt = (s- 1)2 k „ , Wr 1 \t\o. ( g — !) 2 


(s + t- 1 ) k Vt 


(1 + u 2 ) k 2 2 fc -!((fc- l )!) 2 


for integers fc > 1 it follows that for s £ U we have 


v 3 ( S ) = (-i) E ^y 


?:=1 


v^F(2fc- 2)! 
2 2 fe ^ 1 ((fc- 1 )!) 


-afc(s — l ) 2 + R(s) 
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where R is another holomorphic function. Taking the Laurent series of T>\ around 
its pole at s = 1 it follows that T >2 = 2 ?i + 22?3 may be written in the form 


p+i 


V 2 (s) = £ 


fk 


9k 


- +r(s) 


-(s-l)^ (s-l)*-* 

where fk, gk and r are all holomorphic in the region 5ft(s) > 1 , and / p +i(l) = g( 1 ). 
Applying Theorem [2] to V 2 we obtain 


1 


lim 

n-¥o o n z {log n) p 


£ ((C(u/v) — /rlogi;) p + (logi; 


2 p «?(l) 


p! 


and by subtracting the corresponding expression for the partial sums of T>i 
(15) " 1 


lim 

n—too Ti 


- -- £ (logv) p ' 2 (C(u/v)-p\ogv) p = 0. 

Me n„ 


We next show that (fl5l) remains valid if the factor (logr >) p / 2 is removed from inside 
the summation and a factor (log n) p ^ 2 is removed from the denominator outside 
the summation. Following [12], for each t > 0 let us define ip(t) := (logf ) -p / 2 and 
n (*) := E(«,„)efi LtJ {logv) p/2 {C{u/v) - p\ogv) p . Using (HU) we may find r > 2 
such that |II(t)| < 2f 2 (logt) p /p for all t > t and hence for all n > r 


II (t)<p'(t)dt 


< 


n (t)ip'(t)dt 




f n 

<K + J t(fogt)^~ l dtt < K-\-n 2 {\ogn)^~ l , 


say. It follows that 
(16) 


lim 


II (i)ip'(t)dt = 0. 


n->oo ?r 2 (log n) p / 2 

Noting that (log v) p ^ 2 {C{u/v) — n logu) = 0 when v = 1, for each n > 3 we obtain 


f n / / 7 / \ \ P f 71 

/ n (t)ip'(t)dt= £ (logv)S [c (-J - /uloguj / <p'(t)dt 

^ (u,v)£Q, n — 1 v 

£ (log v)% (c 0 ) - £tlog vY ((log n)~% - (logu)" 


(u,v)e o„_i 

1 

(log n) p ! 2 


II(n — 1) — £ (C(u/v) — iilogi’Y 

(u,v)e n„_i 


and in combination with (1151) and (1161) it is not difficult to see that this implies 

1 


(17) 


lim - 

n ->00 n 2 (logn ) p / 2 


£ (C(u/v) — plogv) p = 0. 

(u,i>)EQ n 


To complete the proof of the theorem we must show that for all p > 1, 


( 18 ) 


1 


lim 

n ->00 zfjzii r 


E 

(u,v)(zQ n 


C(u/v) — /ilogn 
Vlog n 


0 if p is odd 

a p (p — 1 )!! if p is even, 
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which differs from (THU) and (fT71) in that the term y log v is replaced with p log n. To 
simplify the notation in the sequel let us write := Er u v )en n (C(u/v) — p log v) p 
and b n , P := E(u,„) e n n (C( u A0 ~ ^ogn) p for each n,p > 1. 

In the case where p is even, let us choose K > 0 such that {plogx) p < K p y/x 
for all x > 1. Using the reverse triangle inequality for the p-norm we obtain 
(19) 





p p jlogrc — log v| 


< K 



Applying Theorem [2] to the series D 0 (s + \) = £“ =1 n 2s (E( M>w ) G n„\n„_ 1 v 1/2 ) 
yields the asymptotic 


( 20 ) 


lim —— 

n—>oo zfjz 


Y J*=(yp) ( lim -±= 4=) 

' V v V fo / \ n ^°° ri\fn ' -Jv I 

(u,v)en n v v ,u 7 \ v v ) 


and in view of (fldll it follows that 

1 

bn.v 


( 21 ) 


lim 


= lim 


‘ri.p 


= o-((p- l)!!) p 


n->oo \#fl n (logn ') p / 2 J n-y oo #fl n (log rc) p / 2 ; 

which yields (1181) for even p. 

We now consider odd p. The difference |a„ jP — b ntP \ may be bounded by 

(22) PH Y, lo gQ (( C Q -^ V Y + ( C (^) _ ^ logn ) P ) 

(n,i>)Ef2 n 


using the elementary inequality \x p — y p \ < p\x— y|(|a;| p_1 + |i/| p_1 ) for real x, y. Let 
e n := (log n)~ 2 / 3 for each n > 1 and consider separately the summation in E^l) over 
pairs ( u , v) such that v < n 1_e ” and over pairs ( u , u) such that v > n 1 ~ Sn . The first 
sum is simply a sum over (u, v) € fii n i- Sn j and so using the bound log(n/u) < log n 
together with m we obtain 

PM (log n) (a|_ n i- e „ JiP _i + b ln i-e n ^ p _ 1 ) = O ((log n) 2 ^ 1 ?! 2-26 ”) = o(n 2 ). 

In evaluating the second sum we may use the estimate log (n/v) < e n log n = 
(logn) 1 / 3 to obtain the upper bound 

pp(logn) 1/3 (a niP _i + 6 n ,p-i) = 0(n 2 (logn)? - £) = o{n 2 (log n) p/2 ). 

We conclude that |a„ iP — b HyP \ = o(n 2 (log?r) p / 2 ) which allows us to deduce m, 
completing the proof for the distributions on Q n . 

In the case of the distribution on , pointwise multiplication of Dirichlet series 
yields 


C(2s)£> p (s) = Y (C(u/v) - plogvf 


n . 

n— 1 d\n \(u,v)eQd\Qd-i 


oo 1 

= E 5 E 


(u,v)^Q n \Q n 


C(u/v) - p log 


gcd(u, v) 
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for 9?(s) > 1, where £ is the Riemann zeta function. To treat this case we consider 
modified Dirichlet series T >\, T > 2,£>3 which are derived from the corresponding series 
T>i byreplacing (C(u/v)—fi log v) with (C(u/v)—/j,log(v/gcd(u,v))) and leaving the 
terms (log-u) p , (logv) p / 2 unchanged. These alternative Dirichlet series are analysed 
by substituting £(2 s)D p (s) for D p (s) throughout the preceding arguments in the 
obvious fashion until we obtain 

i im ± y ( cW*0-/*iog(i5fa ) \ P = ( 0 if P is odd 

\ ^J\ogn I I cr p (p- 1)!! if pis even. 


To derive the analogue of (1181) for even p we follow the model of (1191) . noting that 
the Dirichlet convolution 


C(2s)A) 



£< 


«- 2s £ 
n=1 (u,u)en„\n„_i 


gcd(M,u) 


yields the estimate 

£ 

(u,v)£Q n 

in lieu of (1201) , which suffices to complete the argument for even p. The argument 
for odd p proceeds unchanged. 
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